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Abstract 

The perturbation theory for calculation of the effective conductiv- 
ity of the plane consisting of pieces of different conductivities is con- 
structed and the convenient diagram technique for this perturbation 
theory is elaborated. It is shown that for the chessboard perturbative 
calculations give results which are in agreement with the well-known 
formula (Jeff = ^a\a2- The components of the tensor of effective con- 
ductivity for the anisotropic three-color chessboard are calculated. It 
is shown that the isotropic (symmetric) part of effective conductivity 
calculated up to the sixth order of perturbation theory satisfies the 
Bruggeman effective medium equation for symmetric three-color struc- 
tures with equally partitioned components. We also consider isotropic 
three-color chessboard with nonequal weights of colors. In this case 
the perturbation theory already in fourth order contradicts the results 
following from the Bruggeman equation for nonequal weights. 
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I. 



INTRODUCTION 



The problem of calculation of the effective conductivity of composite mathe- 
rials is of interest from theoretical and phenomenological points of view and 
has been attracting the attention of theorists beginning in the nineteenth' 
century The most developed is the theory of the effective conductiv- 
ity of the plane, because here one can use the property of duality between 
conductivity and resistance, which is typical only for the two-dimensional 
Ohm's law. This duality consists in the fact that the rotation of the plane on 
7r/2 simultaneously with interchange between the vectors of current density 
j and the electric field E does not change Ohm's law. The most interesting 
two-dimensional result is that by Dykhne 0, who has considered the plane 
covered by regions with two different conductivities ai and a2- If the dis- 
tribution is stochastical and the statistical weights of two conductivities are 
equal then one can show that 



The same result will hold also for some regular two-conductivities (for the 
sake of convenience we shall call them below two-color) structures on the 
plane, in particular, for the chessboard as was shown by Keller in his earlier 
paper §. 

The main goal of this paper is the construction of the perturbation theory 
for calculation of the effective conductivity on the plane for an arbitrary dis- 
tribution of the conductivity and the representation of the obtained formulae 
in a convenient graphical form. As applications of this theory we rederive the 
formula of Keller-Dykhne for the case of two-color chessboard. Then, using 
our diagram technique we calculate the components of the tensor of effective 
conductivity for the three-color chessboard which represents an anisotropic 
structure with equal weights. It is shown that the isotropic part of the ef- 
fective conductivity calculated up to the sixth order of perturbation theory 
satisfies the Bruggeman effective medium equation for symmetric three-color 
structures with equally partitioned components ^. 

We consider also the isotropic three-color chessboard with nonequally 
partitioned components. It is shown that in this case the corresponding 
Bruggeman equation failed. 

Notice, that recently, numerical simulation of the conductivity of regular, 
isotropic equal-weighted three-color structures on the plane was performed 



(1.1) 
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and analyzed in the paper 0. Three-component dielectric media were con- 
sidered also in p. 

The structure of the paper is the following: in the second section we 
shall construct the perturbation theory for the effective conductivity and 
present the convenient diagrammatic representation for this technique; the 
third section is devoted to the application of this technique to the deduction of 
the formula for the effective conductivity of the two-color chessboard; finally, 
in Sec. IV we apply perturbative technique for the calculation of components 
of the tensor of effective conductivity for the three-color two-dimensional 
chessboard and show that its isotropic part satisfies the Bruggeman equation, 
then we consider isotropic three-color chessboard with nonequal weights and 
show that fourth-order calculations for this case contradict the predictions 
following from the corresponding Bruggeman equation. 

II. PERTURBATION THEORY FOR THE 
EFFECTIVE CONDUCTIVITY 

Now, let us begin constucting the perturbation theory for the effective con- 
ductivity. The idea of using the perturbative methods for the calculation 
of effective conductivity was considered before by C. Herring [|^ and also 
Dykhne |§] and Bergman 0. Notice that Herring has obtained the formula 
for the tensor of effective conductivity of a locally anisotropic medium in the 
second order 0, while considering locally isotropic mediums we develop the 
technique applicable to arbitrary orders of perturbation theory. 

Ohm's law for a locally isotropic medium has the following form: 

j = aVip, (2.1) 

where j is the current density, a is the conductivity, which generally depends 
on coordinates [x, y) and cp is an electric potential. The charge conservation 
rule for a stationary distribution of currents reads 

divj = V-/=0. (2.2) 

Substituting Eq. ( pTT]) into Eq. ( pl2D one gets 

Av9 + Vlna- Vv3 = 0. (2.3) 
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Let us suppose now that 



a{x,y) = 1 + a{x,y), 
{a) = 0. 



(2.4) 



Here, the average value of a{x, y) is choosen equal to 1 for the sake of conve- 
nience and a{x,y) is a small function whose average on the plane, denoted 
by () is equal to zero. 

Let us represent the potential in the form 



if = E ■ r + ip{x, y). 



(2.5) 



Here, the vector field E should be choosen in such a way to provide the 
fulfillment of the relation 

iyi)) = 0. (2.6) 

Now for isotropic distributions of conductivity one can introduce the effective 
conductivity by means of relation 



where J = (j). For anisotropic distributions of a{x,y) the relation (|2 
becomes a tensor one: 



(2.7) 



Ji = (^effijEj, i,j = 1,2. 



(2.8) 



Meanwhile, substituting the expression for (p from ( |2.5| ) into Eq. ( |2.3| ) 
one gets Laplace equation 



+ V ln(l + a)-{E + Vip) = 0. 



(2.9) 



Now it is convenient to introduce Fourier expansions for the small func- 
tions a and ip: 

a = ^ a^exp(iA; ■ r), 

k 



r . 



(2.10) 
(2.11) 



Substituting expansions ( |2.10D -( p.llD into Eq. ( |2.9D one has 

„2 



V a 



a 



4 
4 



+ 



{E + W) 



(2.12) 



Resolving Eq. ( p.l2D by iterations one can find a perturbative expansion for 
tlie function ipj:. Then, substititing the formulas found for ip^ into Eqs. (|2.7| ) 
or (|2.8|) one can find the components of the tensor of effective conductivity. 



Let us forget for a while about the condition (|2.6|) , providing the correct 
definition of the macroscopic vector E. One can write down the following 
formulas for the perturbative corrections to potenial ipj: and for the isotropic 
part of the tensor of effective conductivity 

(^eff = ]^{(yxx+ (yyy) (2.13) 

i',. = (2.14) 
ik ■ 

^e//2 = - E (2.15) 
k 

^ {kE){pE){pk) 

<^effZ = 1^ ^2^-2^2 "P^g-p^-g' (2-17) 
p,k " 

and so on. 

Making induction one can show that 

, T^n+l^r r {klE){kik2) ■ ■ ■ {kn-2kn^l){kn-lE) 
^e//n - (-1) 2^ fcl, ■ ■ ■ , kn-l ■ ■ ■ k^^, 

Thus all the corrections to effective conductivity have a "chain form" and 
can be represented graphically as 

E ki k2 kn^2 kn-l E 



However, to provide the condition (|2.6|) it is necessary to subtract from 
this diagram all the disconnected diagrams of nth order, consisting of 2 
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chains, then to add the diagrams, consisting of 3 chains and so on. The 
construction of such an alternating sum corresponds to ehmination of the 
terms containing the corrections, proportional to possible homogeneous vec- 
tor harmonics Vip. 

Now, we have to analyze in some detail the structure of expression ( p.l^ ) 
and the corresponding diagram. This diagram contains "vertexes", corre- 
sponding to the coefficients of the Fourier expansion ( p.lOj ), and "propaga- 
tors", corresponding to two-dimensional wave vectors k. Propagators appear 
in 3 forms: 



2' -2fc2+fc2' 



G 



It is remarkable that the formula ( p.l8|) (suplemented, as was explained 
earlier, by alternatingly subtracted terms) contains information about all the 
elements of the tensor of effective conductivity. Diagrams, corresponding to 
the isotropic part of this tensor given by the definition ( |2.13| ) contain an even 



number of propagators F and an even number of propagators G. Diagrams, 
corresponding to anisotropic part of the diagonal elements of this tensor 

^anysot — '^is^xx ^yy)i (2.19) 

contain an odd number of propagators F and an even number of G. Finally, 
diagrams, corresponding to non-diagonal elements of the effective conductiv- 
ity tensor axy, contain an odd number of propagators G. 

III. PERTURBATIVE CALCULATIONS FOR 
THE TWO-COLOR CHESSBOARD 

Now, to show how this perturbation theory and diagram technique works, 
we shall consider the chessboard case. Let us suppose that 

(Ti = 1 + 5, 
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a2 = I - 6. 



(3.1) 



Then 



a. 



2 8 16 



128 



15 



10 



256 



+ 



(3.2) 



To reproduce this formula by means of perturbation theory, we should 
calculate the Fourier coefficients for the function a{x,y). Let us a{x,y) is a 
periodical function with periods 27r in both coordinates which behaves as 



a{x,y) 
a{x,y) 
a{x,y) 



—6, if < X < IT, < y < 7i; 

+6, if n < X < 27r, < y < n; 

+6, if < X < 71, n < y < 2-n; 

—5, if n < X < 27r, n < y < 2tx. 



(3.3) 



y 

27r 

TT 





Then 





-6 


-6 


6 



TC 



271 X 



(y.2k,m = Ctn,2l = 0, (3.4) 

- (27r)2(2/t+l)(2/+l)- ^^-^^ 

It is easy to see that all the diagrams with an odd number of vertices 
give zero contribution. Indeed, the sum of wave numbers of all the vertices 
of the diagram should be equal to zero (see Eq. ( ^.18D ) and only odd values 
of wave numbers are possible (see Eqs. ( |3.4| )- ([375| )). However, the sum of 
an odd number of odd numbers cannot be equal to zero. Thus, we should 
consider only even orders of the perturbation theory. 

Then, if the diagram under consideration contains an odd number of 
propagators G, the contribution of this diagram is equal to zero, because 
for any set of momenta k = kx,ky there is the set k' = —kx,ky, giving 
the contribution of the opposite sign. Analogously, the diagrams containing 
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an odd number of propagators F give vanishing contributions because, the 
contribution of any diagrams with a certain set of momenta k = kx, ky is 
cancelled by that with momenta k = ky^k^. Thus, (Tanisot and a^y are equal 
to zero and the effective conductivity of the chessboard is isotropic. 

Now, we are in a position to calculate the diagrams corresponding to the 
effective conductivity of the chessboard. In the second order of the pertur- 
bation theory one has only the diagram 



which corresponds to the following expression: 

1 



k 

25g^2 oo oo 



1 



2(27r)4 



(3.6) 



Using the well-known formula [10 



oo 1 2 

6 



(3.7) 



one has 



' 2 ■ 



(3.8) 



The set of diagrams of the fourth order has the following form: 



G 



G 



2 •- 



G 



G 



2 •- 



8 



r4 

Here the first connected diagram gives the contribution -g- Combination 



of the second and third diagrams gives again y, the disconnected diagram 
gives the contribution while the last two connected diagrams do not give 
a contribution. Thus, the general result is equal to 
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(3.9) 



In the following picture we shall present only the diagrams, giving non- 
vanishing contributions in the sixth order of perturbation theory with the 
corresponding results: 



G 



G 



G 



G 



F 



G 



G 



G 



G 



G 



G 



G 



G 



32 



32 



_1_ 
32 



• • • 16 



2 ^ ^ 

^« • • • • • 

r) G G 

^« • • • • • 



16 



TE + 



• • • • 



Making summation, one has: 

16 

Quite analogously, one can get also 



^e//6 = -T^- (3.10) 



56' , , 

reffs = -j^ (3.11) 
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= " 256 • ^^-^^^ 
One can see that the results of perturbative calculations ( |3.8| )- (|3TT^ ) coincide 
with those obtained by expansion of the general formula (|3.2|) . It confirms 
the reliability of the perturbation theory developed here and encourage us to 
apply it to more complicated situations. 



IV. PERTURBATIVE CALCULATIONS FOR 
THE THREE-COLOR CHESSBOARDS 

Now we would like to consider three-color structures with square lattice sym- 
metry, which we call the three-color chessboards. To do it, let us turn again 
to the formula (|1 . Ij ) . It is easy to notice that Eq. ( |1 . 1|) could be rewritten in 
the following form: 

The equation ( [4.1|) was written by Bruggeman as an approximative one 
(in the dipole approximation), but appears to be exact for the large class of 
two-color coverings 0]. One can try to generalize this equation for the case 
of three conductivities (three-color system) with equal weights 

which is equivalent to 

^eff + ^(^1 +^2 + fT3)(Te// 

-^(0^10-2 + 0-1 0^2 + Cr2(^3)(^eff " (^l(^2(^3 = 0. (4.3) 

There are different tesselations of the plane allowing three-color (i.e. 
three-conductivities) covering (see, for example, Ref. |ll|). Here we shall 
consider the simplest tesselation: it is again a chessboard but covered by 
three colors with the same statistical weights. We shall introduce again the 
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function a{x,y) which is periodic with periods 27r and is defined by the fol- 
lowing rules: 

27r 27r 

a{x, y)^5, i/ < X < y , < y < y ; 

, , 2n An ^ 2% 

a{x, y) = 0, z/ Y < x < y , < y < y; 

, , - 47r ^ ^ 27r 

y) = Y < a:^ < 27r, < y < y ; 

(^{X: y)^-0, ^/ < X < y, Y < y < y ; 

, , ^ 27r An 2n An 

a{x, y)=d, tf — <x < <y < 

, , ^ „ 47r ^ 27r An 

a{x, y) = 0, y < X < 27r, y < y < y ; 

27r An 

a{x, y) =0, i/ < a; < y , y < y < 27r; 

y) = -(^, U y < a: < y , y < y < 27r; 

An An 
a{x, y)=6, if — <x <2n,— <y <2n. (4.4) 



y, 

27r 



47r 
3 



27r 
3 










-5 


5 




() 





5 





-5 



27r 47r 

T - 27r 



Apparently, this function corresponds to the periodic distribution of three 
conductivitites 



0-2 
0-3 



1+5, 

1-5, 
1. 



(4.5) 
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Correspondingly, the Fourier coefficients have the following form: 

a3k,m = an,3i = 0, (4.6) 

«3fc+l,3i+l = Oi3k+2,3l+2 = 0, (4.7) 
«3fc+2.3.+l = «^,+l,3fc+2 = (2^)2(3^ + 2)(3/ + l)- ^^-^^ 

Now, we are in a position to calculate perturbative contributions to the 
isotropic part of the effective conductivity for three-color chessboard. In the 
second order of perturbation theory we have again only one diagram, the 
contribution of which is equal to 

c^e//2 = -y. (4.9) 



(Here again the formula (|3.71) was used). In the fourth order we have the 



following set of diagrams which have nonvanishing contributions: 
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G G 
• • • • 

I ■ 18 • • • 

• ^ • • • 

which result in 

^e//4 = (4.10) 

Calculating the sums corresponding to the diagrams presented above we have 
used the formula ||T0[ 

S= V - V ^ - V ^ - ^ 

(4.11) 

Below we present the set of all the sixth-order diagrams giving nonvan- 
ishing contributions: 
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G 




G 
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G 


F 




F 


G 


F 


G 




G 


F 




F 




F 






G 




G 





48 



1_ 
72 



1 
108 



1 
216 



G 



F 



G 

• • • 



36 ^ 



• • • • 



1 

108 



Making summation of these contributions, one has 

1756 



<^e//6 - 



432 



(4.12) 



Thus, we have calculated a symmetric part of the tensor of effective con- 
ductivity (Teff = agym = <^xx = <^yy 

Now we can substitute the expression for the isotropic (symmetric) part 
of the effective conductivity found up to the sixth order of the theory of 
perturbation 

(Jeff = 1 (4.13) 

" 3 12 432 ^ ^ 
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into the three-color Bruggeman equation ( [4 .21 ) and to see that it is satisfied. 



It gives us some grounds to think that this equation describes exactly the 
effective conductivity for the chaotic three-color distributions of conductivity 
on the plane just like it does for the two-color distributions, where the formula 
of Bruggeman ( [4.1D coincides with the exact formula of Dykhne. 

It is necessary to stress that for the chessboard with conductivity de- 



scribed by function a{x,y) from Eq. ([4 .4]) there is an anisotropic term al- 
ready in the second order of perturbation theory, which corresponds to the 
following diagram 

G 



and is equal to 

486^2 ^ 1 „ 

{2n)\^t-^{3k + 2){3l + m3k + 2y + {3l + iy) ' ' 

(4.14) 

We would like to stress that for isotropic structures the conductivity in 
the second order is purely local and has universal form which does not depend 
on its particular structure. In anisotropic structures already in the second 
order we encounter non-local propagators and the results is not universal. 

Now let us consider another three-color chessaboard with nonequal weights, 
which however is an isotropic one in contrast to the preceding one. The con- 
ductivity here is given by the formulae 



a{x, y) = 0, if < X < TV, < y < it; 

a{x, y) = +6, ifTT<x< 2n, < y < n; 

a{x, y) = —6, if < X < TV, n < y < 2n; 

a{x,y) =0, if n < X < 2n, n <y <2n. (4.15) 
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y 

27r 



-5 







TT 



2tt X 



Though this chessboard is a three-color one, it possesses all the symme- 
tries of the two-color chessboard and is an isotropic one. 
Non-zero Fourier coefficients have the following form: 



O2fc+l,0 



i7T{2k + l)' 



(4.16) 



It is convenient to denote these Fourier coefficients by different vertices: the 
coefficient from the formula ( [4.16|) by a black circle and that one from formula 
( [4.17| ) by a white circle. 

The contribution of the second order of the perturbation theory into ef- 
fective conductivity is given by the diagram 



and is equal to 

The diagrams of the fourth order have the following form 
• • • 



' 4'' 



(4.18) 



16 



-e- 



-e- 



32 



15 



Thus, the fourth-order contribution is 



5^ 

^eff, = (4.19) 
The Bruggeman equation for this structure could be written as 

2^£ZI^ + ^^ILZ^ + ^^ILZ^ = 0, (4.20) 

or equivalently 

^e// + ^ ^e// ^ ^e// " (TKTzas = 0. (4.21) 

For our case when, 



ai = 1, a2 = 1 + 5, (73 = 1 - 5, 



one can easily get 



-•// = i-T-I5-"-- ("■22' 

It is obvious that in the fourth order of the perturbation theory there is 
a contradiction between the results of perturbative calculations ([4.19|) and 
those following from the Bruggeman equation ( [4.22D . 



ACKNOWLEDGEMENTS 

This work was partially supported by RFBR via grant No 99-02-18409. Work 
of A.K. was also supported by CARIPLO Scientific Foundation. We are 
grateful to A.M. Dykhne, D.J. Bergman and L.G. Fel who have introduced 
us to the problem of effective conductivity of composite systems. We would 
thank J.B. Keller for the interest in our work. We are grateful also to A.S. 
loselevich, I.V. Kolokolov, V.G. Marikhin, A.V. Marshakov, A.Yu. Morozov, 
A. A. Rosly and A.V. Zabrodin for the useful remarks. 



16 



References 

[I] Lord Rayleigh, Phil. Mag. Ser. 5, 34, 481 (1892). 
[2] A.M. Dykhne, Sov. Phys. JETP 32, 63 (1970). 
[3] J.B. Keller, J. Math. Phys. 5, 548 (1964). 
[4] D.A.G. Bruggeman, Ann. Physik (Leipzig) 24, 636 (1935). 
[5] L.G. Fel, V.Sh. Machavariani and D.J. Bergman, |cond- mat/0007022 
[6] Yu.P. Emets, JETP 87, 612 (1998). 

[7] C. Herring, J. Appl. Phys. 31, 1939 (1960). 
[8] A.M. Dykhne, private communication. 
[9] D.J. Bergman, private communication. 

[10] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions, 
New York, Dover, 1965. 

[II] P. Engel, Geometric crystallography: an axiomatic introduction to crys- 
tallography, Kluwer Academic Publishers, Dordrecht, 1986. 



17 



